Detecting Variability in Massive Astronomical Time-Series Data I: application of an infinite Gaussian mixture model

We present a new framework to detect various types of variable objects within massive astronomical time-series data. Assuming that the dominant population of objects is non-variable, we find outliers from this population by using a non-parametric Bayesian clustering algorithm based on an infinite GaussianMixtureModel (GMM) and the Dirichlet Process. The algorithm extracts information from a given dataset, which is described by six variability indices. The GMM uses those variability indices to recover clusters that are described by six-dimensional multivariate Gaussian distributions, allowing our approach to consider the sampling pattern of time-series data, systematic biases, the number of data points for each light curve, and photometric quality. Using the Northern Sky Variability Survey data, we test our approach and prove that the infinite GMM is useful at detecting variable objects, while providing statistical inference estimation that suppresses false detection. The proposed approach will be effective in the exploration of future surveys such as GAIA, Pan-Starrs, and LSST, which will produce massive time-series data.Comment: accepted for publication in MNRA

Broadband Multi-wavelength Properties of M87 during the 2017 Event Horizon Telescope Campaign

In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass ∼6.5 × 10^{9} Mo. The EHTC also partnered with several international facilities in space and on the ground, to arrange an extensive, quasi-simultaneous multi-wavelength campaign. This Letter presents the results and analysis of this campaign, as well as the multi-wavelength data as a legacy data repository. We captured M87 in a historically low state, and the core flux dominates over HST-1 at high energies, making it possible to combine core flux constraints with the more spatially precise very long baseline interferometry data. We present the most complete simultaneous multi-wavelength spectrum of the active nucleus to date, and discuss the complexity and caveats of combining data from different spatial scales into one broadband spectrum. We apply two heuristic, isotropic leptonic single-zone models to provide insight into the basic source properties, but conclude that a structured jet is necessary to explain M87's spectrum. We can exclude that the simultaneous γ-ray emission is produced via inverse Compton emission in the same region producing the EHT mm-band emission, and further conclude that the γ-rays can only be produced in the inner jets (inward of HST-1) if there are strongly particle-dominated regions. Direct synchrotron emission from accelerated protons and secondaries cannot yet be excluded

The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2

We present and analyze observations of polarized dust emission at 850 μm toward the central 1 × 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields in Star-forming Region Observations survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis–Chandrasekhar–Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from Herschel data and the C18O (J = 3 - 2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 ± 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 ± 0.02. Additionally, the mean Alfvén Mach number is 0.35 ± 0.01. This suggests that, in Mon R2, the magnetic fields provide resistance against large-scale gravitational collapse, and the magnetic pressure exceeds the turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically subcritical

Dust polarized emission observations of NGC 6334: BISTRO reveals the details of the complex but organized magnetic field structure of the high-mass star-forming hub-filament network

Context. Molecular filaments and hubs have received special attention recently thanks to new studies showing their key role in star formation. While the (column) density and velocity structures of both filaments and hubs have been carefully studied, their magnetic field (B-field) properties have yet to be characterized. Consequently, the role of B-fields in the formation and evolution of hub-filament systems is not well constrained. Aims. We aim to understand the role of the B-field and its interplay with turbulence and gravity in the dynamical evolution of the NGC 6334 filament network that harbours cluster-forming hubs and high-mass star formation. Methods. We present new observations of the dust polarized emission at 850 μm toward the 2 pc × 10 pc map of NGC 6334 at a spatial resolution of 0.09 pc obtained with the James Clerk Maxwell Telescope (JCMT) as part of the B-field In STar-forming Region Observations (BISTRO) survey. We study the distribution and dispersion of the polarized intensity (PI), the polarization fraction (PF), and the plane-of-The-sky B-field angle (χB_POS) toward the whole region, along the 10 pc-long ridge and along the sub-filaments connected to the ridge and the hubs. We derived the power spectra of the intensity and χBPOS along the ridge crest and compared them with the results obtained from simulated filaments. Results. The observations span 3 orders of magnitude in Stokes I and PI and 2 orders of magnitude in PF (from 0.2 to 20%). A large scatter in PI and PF is observed for a given value of I. Our analyses show a complex B-field structure when observed over the whole region ( 10 pc); however, at smaller scales (1 pc), χBPOS varies coherently along the crests of the filament network. The observed power spectrum of χBPOS can be well represented with a power law function with a slope of-1.33 ± 0.23, which is 20% shallower than that of I. We find that this result is compatible with the properties of simulated filaments and may indicate the physical processes at play in the formation and evolution of star-forming filaments. Along the sub-filaments, χBPOS rotates frombeing mostly perpendicular or randomly oriented with respect to the crests to mostly parallel as the sub-filaments merge with the ridge and hubs. This variation of the B-field structure along the sub-filaments may be tracing local velocity flows of infalling matter in the ridge and hubs. Our analysis also suggests a variation in the energy balance along the crests of these sub-filaments, from magnetically critical or supercritical at their far ends to magnetically subcritical near the ridge and hubs. We also detect an increase in PF toward the high-column density (NH2 â 1023 cm-2) star cluster-forming hubs. These latter large PF values may be explained by the increase in grain alignment efficiency due to stellar radiation from the newborn stars, combined with an ordered B-field structure. Conclusions. These observational results reveal for the first time the characteristics of the small-scale (down to 0.1 pc) B-field structure of a 10 pc-long hub-filament system. Our analyses show variations in the polarization properties along the sub-filaments that may be tracing the evolution of their physical properties during their interaction with the ridge and hubs. We also detect an impact of feedback from young high-mass stars on the local B-field structure and the polarization properties, which could put constraints on possible models for dust grain alignment and provide important hints as to the interplay between the star formation activity and interstellar B-fields

Monitoring the Morphology of M87* in 2009-2017 with the Event Horizon Telescope

The Event Horizon Telescope (EHT) has recently delivered the first resolved images of M87*, the supermassive black hole in the center of the M87 galaxy. These images were produced using 230 GHz observations performed in 2017 April. Additional observations are required to investigate the persistence of the primary image feature- A ring with azimuthal brightness asymmetry- A nd to quantify the image variability on event horizon scales. To address this need, we analyze M87* data collected with prototype EHT arrays in 2009, 2011, 2012, and 2013. While these observations do not contain enough information to produce images, they are sufficient to constrain simple geometric models. We develop a modeling approach based on the framework utilized for the 2017 EHT data analysis and validate our procedures using synthetic data. Applying the same approach to the observational data sets, we find the M87* morphology in 2009-2017 to be consistent with a persistent asymmetric ring of ∼40 μas diameter. The position angle of the peak intensity varies in time. In particular, we find a significant difference between the position angle measured in 2013 and 2017. These variations are in broad agreement with predictions of a subset of general relativistic magnetohydrodynamic simulations. We show that quantifying the variability across multiple observational epochs has the potential to constrain the physical properties of the source, such as the accretion state or the black hole spin

The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament

We present 850 μm dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the subfilaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1-10 pc scales. The magnetic fields revealed in the Planck data are well-aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6-1.0 mG in the DR21 filament and ∼0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by James Clerk Maxwell Telescope. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and subfilaments

First BISTRO Observations of the Dark Cloud Taurus L1495A-B10: The Role of the Magnetic Field in the Earliest Stages of Low-mass Star Formation

We present BISTRO Survey 850 μm dust emission polarization observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the James Clerk Maxwell Telescope (JCMT). We observe a roughly triangular network of dense filaments. We detect nine of the dense starless cores embedded within these filaments in polarization, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field dominated to being matter dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 yr

The Event Horizon Telescope Image of the Quasar NRAO 530

We report on the observations of the quasar NRAO 530 with the Event Horizon Telescope (EHT) on 2017 April 5−7, when NRAO 530 was used as a calibrator for the EHT observations of Sagittarius A*. At z = 0.902, this is the most distant object imaged by the EHT so far. We reconstruct the first images of the source at 230 GHz, at an unprecedented angular resolution of ∼20 μas, both in total intensity and in linear polarization (LP). We do not detect source variability, allowing us to represent the whole data set with static images. The images reveal a bright feature located on the southern end of the jet, which we associate with the core. The feature is linearly polarized, with a fractional polarization of ∼5%-8%, and it has a substructure consisting of two components. Their observed brightness temperature suggests that the energy density of the jet is dominated by the magnetic field. The jet extends over 60 μas along a position angle ∼ −28°. It includes two features with orthogonal directions of polarization (electric vector position angle), parallel and perpendicular to the jet axis, consistent with a helical structure of the magnetic field in the jet. The outermost feature has a particularly high degree of LP, suggestive of a nearly uniform magnetic field. Future EHT observations will probe the variability of the jet structure on microarcsecond scales, while simultaneous multiwavelength monitoring will provide insight into the high-energy emission origin

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimetre wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to 10–100 gravitational radii (rg ≡ GM/c2) scales in nearby sources1. Centaurus A is the closest radio-loud source to Earth2. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our Galactic Centre. A large southern declination of −43° has, however, prevented VLBI imaging of Centaurus A below a wavelength of 1 cm thus far. Here we show the millimetre VLBI image of the source, which we obtained with the Event Horizon Telescope at 228 GHz. Compared with previous observations3, we image the jet of Centaurus A at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that the source structure of Centaurus A resembles the jet in Messier 87 on ~500 rg scales remarkably well. Furthermore, we identify the location of Centaurus A’s SMBH with respect to its resolved jet core at a wavelength of 1.3 mm and conclude that the source’s event horizon shadow4 should be visible at terahertz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses5,6

The TOP-SCOPE Survey of Planck Galactic Cold Clumps: Survey Overview and Results of an Exemplar Source, PGCC G26.53+0.17

This is the final version. Available from American Astronomical Society via the DOI in this record.The low dust temperatures (<14 K) of Planck Galactic cold clumps (PGCCs) make them ideal targets to probe the initial conditions and very early phase of star formation. "TOP-SCOPE" is a joint survey program targeting ∼2000 PGCCs in J = 1-0 transitions of CO isotopologues and ∼1000 PGCCs in 850 μm continuum emission. The objective of the "TOP-SCOPE" survey and the joint surveys (SMT 10 m, KVN 21 m, and NRO 45 m) is to statistically study the initial conditions occurring during star formation and the evolution of molecular clouds, across a wide range of environments. The observations, data analysis, and example science cases for these surveys are introduced with an exemplar source, PGCC G26.53+0.17 (G26), which is a filamentary infrared dark cloud (IRDC). The total mass, length, and mean line mass (M/L) of the G26 filament are ∼6200 M, ∼12 pc, and ∼500 Mpc-1, respectively. Ten massive clumps, including eight starless ones, are found along the filament. The most massive clump as a whole may still be in global collapse, while its denser part seems to be undergoing expansion owing to outflow feedback. The fragmentation in the G26 filament from cloud scale to clump scale is in agreement with gravitational fragmentation of an isothermal, nonmagnetized, and turbulent supported cylinder. A bimodal behavior in dust emissivity spectral index (β) distribution is found in G26, suggesting grain growth along the filament. The G26 filament may be formed owing to large-scale compression flows evidenced by the temperature and velocity gradients across its natal cloud.German Research FoundationJoint Research Fund in AstronomyTop Talents Program of Yunnan ProvinceAcademy of FinlandMinistry of Education, Science, and TechnologyNational Research Foundation of KoreaChinese Academy of SciencesMinistry of Science and Technology of TaiwanEuropean Research Counci